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Abstract 

In this paper, the effect of electromagnetic field has been investi- 
gated on the spherically symmetric collapse with the perfect fluid in 
the presence of positive cosmological constant. Junction conditions 
between the static exterior and non-static interior spherically sym- 
metric spacetimes are discussed. We study the apparent horizons and 
their physical significance. It is found that electromagnetic field re- 
duces the bound of cosmological constant by reducing the pressure 
and hence collapsing process is faster as compared to the perfect fluid 
case. This work gives the generalization of the perfect fluid case to the 
charged perfect fluid. Results for the perfect fluid case are recovered. 
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1 Introduction 

General Relativity (GR) predicts that gravitational collapse of massive ob- 
jects (having mass = IO^Mq — 10^ Mq, where Mq is mass of the Sun) results 
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to the formation of spacetime singularities in our universe pQ . The singularity 
theorems [2] of Hawking and Penrose reveal that if a trapped surface forms 
during the collapse of compact object, such a collapse will develop a space- 
time singularity. According to these theorems, the occurrence of spacetime 
singularity (can be observed or not) is the generic property of the spacetime 
in GR. An observable singularity is called naked singularity while other is 
called black hole or covered singularity. 

An open and un-resolved problem in GR is to determine the final fate of 
the gravitational collapse (i.e., end product of collapse is either covered or 
naked singularity). To resolve this problem, Penrose p] suggested a hypoth- 
esis so-called Cosmic Censorship Conjecture (CCC). This conjecture states 
that the singularities that appear in the gravitational collapse are always 
covered by the event horizon. It has two versions, i.e., weak and strong 
version [1]. The weak version states that the gravitational collapse from 
the regular initial conditions never creates spacetime singularity visible to 
distant observer. On the other hand, the strong version says that no singu- 
larity is visible to any observer at all, even some one close to it. There is no 
mathematical or theoretical proof for either of the version of the CCC. 

The singularity (at the end stage of the gravitational collapse) can be 
black hole or naked depending upon the initial data and equation of the 
state. To prove or disprove this hypothesis, many efforts have been made 
but no final conclusion is drawn. It would be easier to find the counter 
example that would enable us to claim that the hypothesis is not correct. 
For this purpose, Virbhadra et al. [5] introduced a new theoretical tool using 
the gravitational lensing phenomena. Also, Virbhadra and Ellis [6| studied 
the Schwarzschild black hole lensing and found that the relativistic images 
would confirm the Schwarzschild geometry close to the event horizon. The 
same authors [7j analyzed the gravitational lensing by a naked singularity 
and classified it into two kinds: weak naked singularity (those contained 
within at least one photon sphere) and strong naked singularity (those not 
contained within any photon sphere). 

Claudel et al. [S] showed that spherically symmetric black holes, with 
reasonable energy conditions, are always covered inside at least one photon 
sphere. Virbhadra and Keeton [9] studied the time delay and magnification 
centroid due to gravitational lensing by black hole and naked singularity. It 
was found that weak CCC can be tested observationally without any ambi- 
guity. Virbhadra [10] explored the useful results to investigate the Seifert's 
conjecture for naked singularity. He found that naked singularity forming in 
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the Vaidya null dust collapse supports the Seifert's conjecture. In a recent 
paper [TT], the same author used the gravitational lensing phenomena to find 
the improved form of the CCC This work is a source of inspiration for many 
leading researchers. 

Oppenheimer and Snyder |12] are the pioneers who investigated gravi- 
tational collapse long time ago in 1939. They studied the dust collapse by 
taking the static Schwarzschild spacetime as exterior and Friedmann like so- 
lution as interior spacetime. They found black hole as end product of the 
gravitational collapse. To study the gravitational collapse, exact solutions of 
the Einstein field equations with dust provide non-trivial examples of naked 
singularity formation. Since the effects of pressure cannot be neglected in 
the singularity formation, therefore dust is not assumed to be a good matter. 

There has been a growing interest to study the gravitational collapse in 
the presence of perfect fiuid and other general physical form of the fiuid. Mis- 
ner and Sharp [13] extended the pioneer work for the perfect fluid. Vaidya 
[H] and Santos [15] used the idea of outgoing radiation of the collapsing body 
and also included the dissipation in the source by allowing the radial heat 
flow. The cosmological constant A affects the properties of spacetime as it 
appears in the fleld equations. It is worthwhile to solve the fleld equations 
with non-zero cosmological constant for analyzing the gravitational collapse. 
Markovic and Shapiro [I6] generalized the pioneer work with positive cos- 
mological constant. Lake P3] extended it for both positive and negative 
cosmological constant. 

Sharif and Ahmad [I8]-[2T] extended the spherically symmetric gravita- 
tional collapse with positive cosmological constant for perfect fluid. They 
discussed the junction conditions, apparent horizons and their physical sig- 
niflcance. It is concluded that apparent horizon forms earlier than singular- 
ity and positive cosmological constant slows down the collapse. The same 
authors also investigated the plane symmetric gravitational collapse using 
junction conditions [22]. In a recent paper [23], Sharif and Iqbal extended 
plane symmetric gravitational collapse to spherically symmetric case. 

Although a lot of work has been done for dust and perfect fluid collapse 
of spherically symmetric models. However, no such attempt has been made 
by including the electromagnetic fleld. We would hke to study the gravita- 
tional collapse of charged perfect fluid in the presence of positive cosmological 
constant. For this purpose, we discuss the junction conditions between the 
non-static interior and static exterior spherically symmetric spacetimes. The 
main objectives of this work are the following: 
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• To study the effects of electromagnetic field on the rate of collapse. 

• To see whether or not CCC is valid in this framework. 

The plan of the paper is as follows: In the next section, the junction 
conditions are given. We discuss the solution of the Einstein-Maxwell field 
equations in section 3. The apparent horizons and their physical significance 
are presented in section 4. We conclude our discussion in the last section. 

We use the geometrized units (i.e., the gravitational constant G=l and 
speed of light in vacuum c=l so that M = ^ and k = = Svr). All the 
Latin and Greek indices vary from to 3, otherwise it will be mentioned. 



2 Junction Conditions 

We consider a timelike 3D hypersurface E which separates two 4D manifolds 
M~ and M"*" respectively. For the interior manifold M~ , we take spherically 
symmetric spacetime given by 

dsl = - X^dr^ - Y^de^ + sin^'d^'), (2.1) 

where X and Y are functions of t and r. For the exterior manifold we 
take Reissner-Nordstrom de-Sitter spacetime 

dsl = NdT^ - ^dR^ - R\d9^ + sin ^2^02), (2.2) 

where 

2M g2 ^ 

NiR) = l- — + ^--R\ (2.3) 

M and A are constants and Q is the charge. The Israel junction conditions 
are the following [24j : 

1 . The continuity of line element over E gives 

{dslh = {dslh = dsl. (2.4) 

2. The continuity of extrinsic curvature over E yields 

[i^.,] = - = 0, (z,j = 0,2,3) (2.5) 
where Kij is the extrinsic curvature defined as 
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^^ = -^^(a^ + r^^^|i^)' (o^,/i,^ = 0,l,2,3). (2.6) 

Here ^* correspond to the coordinates on S, stand for coordinates in 
M^, the Christoffel symbols are calculated from the interior or exterior 



spacetimes and are components of outward unit normals to S in the 
coordinates xZ^. 

The equation of hypersurface in terms of interior spacetime M' coordi- 
nates is 

/-(r,t) = r-rs = 0, (2.7) 

where is a constant. Also, the equation of hypersurface in terms of exterior 
spacetime coordinates is given by 

/+(/?, r) = i?-i?s(T) = 0. (2.8) 

When we make use of Eq. (12.71) in Eq. (l2.1l) . the metric on S takes the form 

{ds'_)Y. = dt^ -Y'^{r^,t){de'^ ^^\^Q'^d(\?). (2.9) 

Also, Eqs.dMl) and ([2j) yield 

{ds\)^ = mn^) - j^^(^-^f]dT' - RUd9' + sme'd<t>'), (2.10) 
where we assume that 

SO that T is a timelike coordinate. From Eqs. fl2.4p . (12.91) and (I2.10p . it follows 

that 

/iE = K(rE,i), (2.12) 

Also, from Eqs. (12.71) and (12. 8p . the outward unit normals in M" and M"*", 
respectively, are given by 

n- = (0,X(r2,t),0,0), (2.14) 
n+ = (-i?s,t,0,0). (2.15) 



The components of extrinsic curvature become 

^00 = 0, (2.16) 

YY' 

K^, = csc'9K^, = {^h, (2.17) 

, NdN 3 dN ■ -y, 

AV„ = (RT-TR---T^^--TR%, (2.18) 

A'i = CSC- eKi, = (NRt)s, (2,19) 

where dot and prime mean differentiation with respect to t and r respectively. 
From Eq.f l2.5p . the continuity of extrinsic curvature gives 

Ko+ = 0, (2.20) 

K+ = K-,. (2.21) 

Using Eqs. fl216D - fl2:2TD along with Eqs.([23D, fl232D and fl27[3D . the junction 
conditions become 



M={l-^y' + ^. + ly"- 7^Y"h- (2-23) 



{XY' - XY% = 0, (2.22) 
2 6 2y 2 2X2 

3 Solution of the Einstein Field Equations 

The Einstein field equations with cosmological constant are given by 

G,,-Ag,, = K{T,, + TlfJ^^). (3.1) 
The energy-momentum tensor for perfect fluid is 

T^,u = (p + p)u^Uu - pg^u, (3.2) 

where p is the energy density, p is the pressure and = S'^ is the four- vector 

velocity in co-moving coordinates. T^l^"*^ is the energy-momentum tensor for 
the electromagnetic field defined [25] as 

Tjr^ = ^{-g'^^F^sF,^ + ^g,,Fs^F'n. (3.3) 
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With the help of Eqs.([S2D and Eq.(^ takes the form 

R,u = 87r[(p + p)u,u, + ^{p- p)g^, + T^^) - \g,.T^'"'^ - Ag,.. (3.4) 

To solve this equation, we need to calculate the non-zero components and 
trace free form of T^t'^\ For this purpose, we first solve the Maxwell's field 
equations 

F^^u = <Pu,f, - (pfi,u, (3.5) 
F'"';iy = -471 J'^, (3.6) 

where is the four potential and J'^ is the four current. As the charged 
fluid is in co-moving coordinate system, the magnetic field will be zero in 
this case. Thus we can choose the four potential and four current as follows 

= (0(t,r), 0,0,0), (3.7) 
J" = au", (3.8) 

where a is charge density. 

Now for the solution of the Maxwell's field Eq. (13.61) . and u are treated 
as local coordinates. Using Eqs. (13.51) and (13. 7p . the non-zero components of 
the field tensor are given as follows: 

Foi = -Fio = (3.9) 



Also, from Eqs. (l3.6p and (13.81) . we have 

1 d'^cp _ d(j) X' 
X dr'^ dr X"^ 

1 d^d) X 



-4:%aX, (3.10) 



^lf = 0- (3-11) 



X drdt X^ dr 
Equation (13.111) implies that 

(If)^*-, (3.12, 

where K = K[r) is an arbitrary function of r. Equations (13.101) and (13.121) 
yield 

K'{r) = -AnaX. (3.13) 
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The non-zero components of T^l^™'-' and its trace free form turn out to be 

rp{em) _ 1 T^2 rp{em) _ 1 t^2 y-2 rpiem) _ 1 T^'2^2 
loo - ' - ^ ' ^22 ^ ' 

r3(|™) = rir^sin^^, t(^-) = o. 

When we use these values, the field equations (13.41) for the interior space- 
time takes the form 

Roo = -^-2y = An{p + 3p) + K^ -A, (3.14) 

X XY 2 Y" Y'X' 
^11 - 'x' XY^X^^Y~1(Y^ 

= A7r{p - p) + - A, (3.15) 

r> ^ \2 XY 2 fY" 1^' \2 X' Y' X 2i 

-fi22 — —y ~ ly/ ~ X y ~'~ 'X^ ~Y' ~'~ r ~ 'Jc~Y ~ r 

= 47r(p - p) - - A, (3.16) 

i?33 = sin2^/222, (3.17) 
Y' X Y' 

Roi = -^y + ^xy = ^- (3-18) 

Now we solve Eqs. fl3.14p - fl3.18p . Integration of Eq. fl3.18p with respect to t 
yields 

X = ^, (3.19) 

where H = H{r) is an arbitrary function of r. The energy conservation 
equation 

T;^, = 0, (3.20) 

for the perfect fluid with the interior metric shows that pressure is a function 
of t only, i.e., 

p = pit). (3.21) 

Substituting the values of X and p from Eqs. fl3.19p and fl3.2ip in Eqs. (13.141) - 
f l3A8|) . it follows that 

Y Y (1 — R'^\ 
2^ + (^)^ + = A + K^- Pit). (3.22) 
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We consider p as a polynomial in t as given by 

Pit) = Pet-', (3.23) 

where Pc and s are positive constants. Further, for simplicity, we take s = 
so that 

Pit) = p,. (3.24) 
Replacing this value in Eq. fl3.22p . we get 

2y + iy) + y2 = A + - 8m. (3.25) 
Integrating this equation with respect to t, it follows that 

Y'^ 7T) 

Y^ = H'-l + (A + K'- 8np,)— + 2y, (3.26) 

where m = m(r) is an arbitrary function of r and is related to the mass of 
the collapsing system. Substituting Eqs. fl3.19p . (13.260 into Eq. fl3.14l) . we get 

9K'K 

m = -^Y^ + y'y'[47r(pe + p) + (3.27) 

For physical reasons, we assume that pressure and density are strictly 
positive. Integrating Eq.f l3.27p with respect to r, we obtain 

m(r) = 47r / (p + Pe)Y'Y^dr + 2 K'^Y'Y^dr + - / K'KY^dr. (3.28) 
Jo Jo 3 Jq 

The function m(r) must be positive because m(r) < implies negative mass 
which is not physical. Using Eqs. fl3.19p and fl3.26p into the junction condition 
Eq.fESS]), it follows that 

O"^ 1 

M = ^ + m + -iA + - 87ipc)Y^ . (3.29) 

The total energy M(r, t) up to a radius r at time t inside the hypersurface 
S can be evaluated by using the definition of mass function [13] given by 

Mir,t) = ^Yil+g^''Y^Y,). (3.30) 
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For the interior metric, it takes the form 



M{r^) = \Y{l + Y^-{^f). (3.31) 
Replacing Eqs. fl3.19l) and fl3.26p in Eg. (13.311) . we obtain 

Mir,t) =m{r) + {K + K'^ -Mp^) — . (3.32) 

6 

Now we take (A + K'^ — SttjOc) > and the assumption 

H{r) = 1. (3.33) 

In order to obtain the analytic solutions in closed form, we use Eqs. (l3.19|l . 
(K26\i and (K33\i so that 

y = ( A + ,f2_8.,/ ^""'«('-.0. (3.34) 

/ 6m ,1 ,^m' 2KK' . . , / x 
X = ^ )^ { 7-. n -}smha(r,t) 

^ 2{ts{r)-t)KK' , , / A + K2-87rpe , 
+ V3(A + i^^-8m)+ ^^'^V 3 ^ 

X cosh a(r, if:)] sinh~ a(r, t), (3.35) 

where 

a(r,t) = ^^ ^ [ts{,r)-t)]. (3.36) 

Here ts(T') is an arbitrary function of r and is related to the time of formation 
of singularity of a particular shell at coordinate distance r. 

In the limit (Syrpc — K'^) A, the above solution corresponds to the 
Tolman-Bondi solution [26j 

hm Xir,t) = ^'(ts-t) + 2rnt' 

hm Y{r,t) = [^(t.-t)^]i (3.38) 
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4 Apparent Horizons 

Here we discuss the apparent horizons for the interior spacetime. The bound- 
ary of two trapped spheres whose outward normals are null is used to find 
the apparent horizons. This is given as follows: 

g^''Y,,Y,, = Y'- (^f = 0. (4.1) 
Replacing Eqs. fl3.19l) and fl3.26p in this equation, we get 

{A + K^-Snp^) 3r + 6m = 0. (4.2) 

3 

When we take A = Snpc — K^, it gives Y = 2m. This is called Schwarzschild 

3 



horizon. For m = pc = K = 0, we have Y = y j^, which is called de-Sitter 
horizon. Equation (14.21) can have the following positive roots. 

Case (i): For 3m < , , we obtain two horizons 

Yi = ^ COS-, (4.3) 

^(A + _ g^pj 3 

Y, = ~\ (cos|- V3sin^), (4.4) 

^y{A + SnK^ - Pc) 3 3 

where 



cos if = -3mv/(A + i^2 _87rj9,). (4.5) 



If we take m = 0, it follows from Eqs. (14.31) and (14. 4p that Yi — y (^x+k'^-%ttp ) 
and 12 = 0. Yi and Y2 are called cosmo logical horizon and black hole horizon 
respectively. For m 7^ and A 7^ Svrpc — K^-, Yi and Y2 can be generalized 
respectively. 



Case (ii): For 3m = , ^ , there is only one positive root which 
corresponds to a single horizon i.e., 

Yi = Y2= ^ = Y. (4.6) 

V(A + K^- 8np,) 
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This shows that both horizons coincide. The range for the cosmological and 
black hole horizon can be written as follows 



The black hole horizon has its largest proper area 

cosmological horizon has its area between ^^_^^2^g^^ ^ and ^^^^^g^^ ^ . 

Case (ill): For 3m > . , there are no positive roots and con- 

sequently there are no apparent horizons. 

We now calculate the time of formation for the apparent horizon using 
Eqs.( l333|) . (Km and M 

ts- , ^ sinh-^(— (n = l,2). (4.8) 



When Sirpc — K — ^ A, this corresponds to Tolman-Bondi [26] 



From Eq. fl4.8p . we can write 



4 

tah = ts- -m. (4.9) 



^ = cosh2a,, (4.10) 
2m 



where a„(r,t) = V3(^+^^ ^"^^^ ■[t,(r) - Q]. Equations g2D and gSD imply 
that Yi > Y2 and ^2 > respectively. The inequality ^2 > indicates that 
the cosmological horizon forms earlier than the black hole horizon. 

The time difference between the formation of cosmological horizon and 
singularity and the formation of black hole horizon and singularity respec- 
tively can be found as follows. Using Eqs. fl4.3l) - fl4.5p . it follows that 

d{^) 1 , sin ^ 3cos^, 

= ^ + ^ <0, 4.11 

dm m sm (p cos (p 

1 sin(^ Scos^^, , , 

= — ( ^ + ^)>0. 4.12 

dm m sin ip cos ip 
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The time difference between the formation of singularity and apparent hori- 
zons is 

Tn = ts-tn. (4.13) 

It follows from Eq. (l4:T0D that 

' (4.14) 



d{^) sinh an cosh a/3 (A + K"^ - 8np^) 
Using Eqs.dlJl]) and dHHD, we get 

dT, dT, di^) 1 



dm d{^) dm m^J?>{K + - Mpc) sinh ai cosh ai 

, sin I 3 cos ^ , , ^ 

X ^ + ^ <0. 4.15 

sm if cos Lp 

It shows that Ti is a decreasing function of mass m. This means that time 
interval between the formation of cosmological horizon and singularity is 
decreased with the increase of mass. Similarly, from Eqs. (l4.12l) and (l4.14p . 

we get 

dT. 1 



dm mY/3(XT~K^"--"87i7)J sinh q;2 cosh 0^2 

sin(^ 3cos^^, , , 

x( ^ + ^)>0. (4.16) 

sin cos Lp 

This indicates that T2 is an increasing function of mass m indicating that 
time difference between the formation of black hole horizon and singularity 
is increased with the increase of mass. 



5 Summary and Conclusion 

This paper is devoted to study the effects of electromagnetic field on gravi- 
tational collapse with the positive cosmological constant. The cosmological 
constant acts as Newtonian potential. The relation for the Newtonian poten- 
tial is = |(1 — (7oo)- Using Eqs. fl2.12l) and (13.291) . for the exterior spacetime, 
the Newtonian potential turns out to be 

</>(/?) = ^ + (A + K2-87rp,):f. (5.1) 
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The corresponding Newtonian force is 

F = -^ + iA + K'-87Tp,)^. (5.2) 

Now we discuss the consequence of the Newtonian force. This force is zero 
for the fixed values of m = — . ^ = and R = , and will be 

3y/{A+K2-87Tpc) i/(A+ii-2-87rpc) 

positive (repulsive) if the values of m and R are taken larger than these values. 
If we take m = , ^ and R = , ^ then F = 2(A+i^^-8^Pc) 

which gives positive value if (A + — Sfrpc) > 0, i.e., A > {Snpc — K"^) such 
that 8npc > . Thus we conclude that the repulsive force can be generated 
from A if A > (Svrpc ~ K"^) such that Syrpc > K"^ over the entire range of 
the collapsing sphere. For the perfect fluid and dust cases, A can play the 
role of the repulsive force for A > Svrpc and A > respectively. Notice that 
K = K{r) gives the electromagnetic field contribution. From Eq. fl3.26p . the 
rate of collapse turns out be 

y = -^ + (A + K2-87rp,)-. (5.3) 

This shows that we have re-formulated the Newtonian model which repre- 
sents the acceleration of the collapsing process. The analysis of positive and 
negative acceleration would give the same results as for the Newtonian force. 

It is worthwhile to mention that the electromagnetic field reduces the 
bound of the positive cosmological constant by reducing the pressure. Thus 
the positive cosmological constant is bounded below as compared to the 
perfect fluid case. This would decrease the repulsive force which slows down 
the collapsing process. Making the analysis of the smaller values of m and 
R than the values used for the repulsive force, we find that the attractive 
force is larger than the perfect fluid case. Since the attractive force favors 
the collapse while the repulsive force resists against the collapse, thus the 
collapsing process is faster as compared to perfect fluid case when we include 
the electromagnetic field. 

Further, we have found two apparent horizons (cosmological and black 
hole horizons) whose area decreases in the presence of electromagnetic field. 
It is found that the cosmological horizon forms earlier than the black hole 
horizon. Also, Eq. fl4.8p shows that apparent horizon forms earlier than sin- 
gularity. In this sense, we can conclude that the end state of gravitational 
collapse is a singularity covered by the apparent horizons (i.e., black hole). 
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It is interesting to mention here that our study supports the CCC and 
would be considered as one of its counter example. Also, it would be possible 
that the electromagnetic field reduces the range of apparent horizons to ex- 
treme limits and singularity would be locally naked. Thus the weak version 
of the CCC seems to be valid in this case. 
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